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ABSTRACT 

We present ten medium-resolution, high signal-to-noise ratio near- infrared (NIR) spectra of 
SN 201 lfe from SpeX on the NASA Infrared Telescope Facility (IRTF) and Gemini Near-Infrared 
Spectrograph (GNIRS) on Gemini North, obtained as part of the Carnegie Supernova Project. This 
data set constitutes the earliest time-series NIR spectroscopy of a Type la supernova (SN la), with the 
first spectrum obtained at 2.58 days past the explosion and covering —14.6 to +17.3 days relative to 
_B-band maximum. C I A1.0693 /im is detected in SN 2011fe with increasing strength up to maximum 
light. The delay in the onset of the NIR C I line demonstrates its potential to be an effective tracer 
of unprocessed material. For the first time in a SN la, the early rapid decline of the Mg II A1.0927 
lim velocity was observed, and the subsequent velocity is remarkably constant. The Mg II velocity 
during this constant phase locates the inner edge of carbon burning and probes the conditions under 
which the transition from deflagration to detonation occurs. We show that the Mg II velocity does not 
correlate with the optical light-curve decline rate /^m\^{B). The prominent break at ~ 1.5 /mi is the 
main source of concern for NIR k-correction calculations. We demonstrate here that the feature has 
a uniform time evolution among SNe la, with the flux ratio across the break strongly correlated with 
Ami 5 (B). The predictability of the strength and the onset of this feature suggests that the associated 
k-correction uncertainties can be minimized with improved spectral templates. 

Subject headings: cosmology: observations — infrared: general — supernovae: general — supernovae: 
individual (SN 201 lfe) 



1. INTRODUCTION 

SNe la provide a direct measure of the expansion 
history of the universe and have l ed to the discov- 
ery of the accelerate d expansion (jRiess et all 119981 : 
iPerlmutter et alJ ll999'). The unknown cause of the ac- 
celerated expansion is commonly referred to as "dark 
energy" . Although SNe la remain the most proven 
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technique for studying dark energy, there are legitimate 
concerns that systematic errors coming from the astro- 
physics of SNe la will ultimately limit their accuracy. 

Fort unately, observatio ns in the NIR offer a way for- 
ward. lElias et all (|1985H presented the first SN la NIR 
Hubble diagram with a dispersion of only 0.2 mag. Sub- 
sequent studies indicated that NIR peak luminosity in- 
deed show a smaller intrinsic scatter than that in the 
optical fe.g.. | Krisciunas et al.l l2004t iWood-Vasey et al.l 
[20081: iFolatelli et al.l 120101: iBarone-Nugent et al.l l2012f) 7 
a result which was supported by the theoretical work 
of iKasenl (|2006f) . As op posed to the width-luminosity 
relation in the optical (jPhillipsI 119931 ). NIR peak lu- 
minosities for all except the fastest-declining events 
have only a slight dependence on the l ight-c urve shape 
(|Krisciunas et al.l 120091: iKattner et ail |2012[ ). Further 
analyses confirmed that NIR luminosities without light- 
curve shape corrections are more effective distance indi- 
cators than corrected optical luminosities (|Mandel et al.l 
120091 120111 ) . A key ingredient to realize the full poten- 
tial of NIR SN la cosmology is NIR spectroscopy, such 
that the peak luminosities can be accurately k-corrected 
to the rest frame. With the limited scope and size of 
the world's current NIR spectroscopic sample, the time 
evolution and the diversity of the spectral features are 
to date poorly understood. These uncertainties directly 
affect the determination of the NIR rest-frame peak lu- 
minosities. 

Given that SNe la play a critical role in observational 
cosmology, it is as critical as ever to understand the ex- 
plosion mechanisms and the progenitor systems. NIR 
spectra also provide several key diagnostics of SN la 
physics. The C I A1.0693 /^m line is strong and iso- 
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lated in the NIR and is an ideal tracer of primordial, un- 
processed material f rom the carbon-oxygen white dwarf 
(|Marion et al.l [20061) . The amount and the distribution 
of unprocessed material are key discriminators for SN la 
explosion models. Magnesium is a product of explosive 
carbon burning and not oxygen burning. The strong NIR 
Mg II A1.0927 (im line provides a direct measure of the 
inner boundary of carbon burning and probes the con- 
ditions under w hich transition from deflagration to det- 
onation occurs (jWheeler et aLlfToM) . After maximum 
light, the dramatic drop in electron scattering opacity, 
coupled with heavy line blanketing from fully-processed 
iron-peak material, rapidly forms the prominent feature 
at 1.5 — 1.9 fim and prov ides a temperature p robe of the 
line blanketing material (jWheeler et al.| [l998). 

With the advent of NIR narrow-gap semiconductor de- 
tectors, pioneering studie s of SN la NIR spectroscopy be- 
gan i n the late 1980s fe.g..lFrogel et al.lll987tlLvnch et all 
119901 ISpvromilio et alj f"l992. 1994ft. IBowers et al.l (119971 ) 
published the first multi-object sample, which consisted 
of NIR spectra taken in late and nebular phases. The cur- 
rent sample of SN la NIR spectroscopy largely consists 
of data of a few nearby SNe la. We list h ere the spec- 
tra in cluded in our analysis: SNe 1994D dMeikle et all 
19961) . 1998bu (iJha et a ll U99l Hernandez et all 120001: 



Hamuv et al.ll2002l). 1 999bv (IHoflich et alll2002h 1999ee 



(lHamuv et al.l 120021) 2002 bo dBenetti et all l200l 
2002d i (IPignata et al.l I2008D. 2003d u (IStanishev et all 
20071), 2004S dKrisciunas et al.ll2007h. 2005cf dGall et all 



2012), 2011iv (|Folev et all 120121 ). iMarion et all (|2003l 
20091 ) published a catalog of 41 SN la near-infrared spec- 



tra, nearly doubling the world's sample. Owing to these 
works, we can begin to examine the statistical spectro- 
scopic properties of SNe la in the NIR. 

Despite the progress made thus far, the NIR spectro- 
scopic sample of SNe la is not yet ideal. The NIR sample 
size is two orders of magnitude smaller than that in the 
optical. Many of the SNe la observed lack accompanying 
photometric information. Most importantly, there has 
only been a handful of time-series observations. Densely- 
sampled time-series observations are required to under- 
stand how the features evolve with time and vary among 
object-to-object. The discovery of the nearby SN 2011fc 
made possible such observations. 

On August 24, 2011, SN 201 lfe was detected withi n 
hours of its explosion in M101 (|Nugent et al.l 1201 lh . 
Its proximity and early detection provided a unique 
opportunity to make exquisitely detailed observations of 
a SN la. SN 20 life has reportedly representative proper- 
ties of SNe la and serves as an ideal baseline to compare 
to other objects. It is rapidly becoming the best-studied 
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explosion images of SN 201 lfe placed strong constraints 
on the progenitor system, ruling out a massive compan- 
ion star, disfavoring most theoretical double-degenerate 
progenitor systems and excluding Roche-Lobe over- 
flowing red giant and main- sequence companions to 
high signific ance (iNu gent et all 120111: iLi et all 120111: 
iBloom et al.l I2012D ; although, a M dwarf companion 



may evade these constraints (Wheeler 20 12}). Ten NIR 
spectra of SN 20 life were obtained in the span of a 
month, as part of a joint Carnegie Supernova Project 
(CSP)-CfA Supernova Group program to obtain a 
statistical sample of NIR spectroscopic observations of 
supernovae (Phillips et al. in prep). The SN 201 lfe data 
set constitutes the earliest and highest signal-to-noise 
ratio time-series NIR spectroscopy of a SN la. In this 
paper, we present several analyses and insights on time 
evolution, which would not have been possible without 
this one-of-a-kind data set. 

2. OBSERVATION AND REDUCTION 

Only 2.58 days past the explosion, the first NIR spec- 
trum of SN 201 lfe was obta ined with the SpeX spec- 
trograph (Ra vner et all[2003l ) on IRTF. Nine subsequent 
NIR spectra were o btained with a re gular three-day ca- 
dence with GNIRS (jElias et al.l 11998 s ) on Gemini North. 
Relatively short per-frame exposure times were chosen 
to prevent the saturation of the telluric OH lines and 
the supernova. The object was nodded along the slit 
using the classical ABBA technique. Two ABBA se- 
quences (8 exposures) were obtained at each epoch. To 
correct the science spectra for the effect of telluric at- 
mospheric absorption, an A0V star was observed at a 
similar air mass to that of the science observation. The 
slit was oriented along th e parallactic angle for all obser- 
vations (|Filippenkolll982l ). The spectra from both SpeX 
and GNIRS presented here have similar medium spectral 
resolving power of R ~ 1000 — 2000. With the exception 
of the first SpeX spectrum, all of the resulting spectra 
have very high signal-to-noise ratio of well above 100. 
A journal of the spectroscopic observations is given in 
Table [U 

The SpeX spectrum was obtained in the cross- 
dispersed mode, utilizing a grating and prism cross- 
dispersers, with a 0'.'5 slit. This configuration yields a 
wavelength coverage from 0.8 to 2.5 /zm divided over six 
orders, with gaps over the regions of strong telluric ab- 
sorptions between the JHK bands. With the choice of 
the 0'.'5 slit width, the resulting spectral resolution is in 
the range of ~ 1600 — 2000. The data were calibrated and 
reduc ed using the IDL pipeline Spextool (Cushing et al. 
I2004} ). specifically designed for the reduction of SpeX 
data. The pipeline performed steps of pair subtraction, 
flat-fielding, aperture definition, spectral tracing and ex- 
traction, and wavelength calibration. The ID spectra 
from separate exposures were then combined. Correc- 
tions for telluric absorption we re performed using the 
IDL tool xtellcor developed bv lVacca et al.l (|2003l) . To 
construct a telluric correction spectrum free of stellar 
absorption features, a model spectrum of Vega was used 
to match and remove the hydrogen lines of Paschen and 
Brackett series from the AV0 telluric standard. The re- 
sulting telluric correction spectrum was also used for the 
absolute flux calibration, as the NIR continua of A0V 
stars are reasonably approximated by a blackbody, with 
the temperature (typically ^10,000 K) determined by 
the color of the star. 

The GNIRS spectra were observed in the cross- 
dispersed mode, in combination with the short- 
wavelength camera, a 32 lines per mm grating and a 
l'.'O slit. This configuration allows for a wide contin- 
uous wavelength coverage from 0.8 to 2.5 /jm, divided 
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TABLE 1 
List of Observations 



UT 


UT 




Phase wrt 


Phase wrt 


Number of 


Integration 


SN 2011fe 


Telluric 


Telluric 


Date 


MJD 


Instrument Explosion 


B Maximum 


Exposures 


Time 


Air mass 


Standard 


Air mass 


2011-08-26 


55799.27 


SpeX 


2.58 d 


-14.6 d 


8 


1200 s 


2.00 


HD143187 


1.92 


2011-08-28 


55801.25 


GNIRS 


4.56 d 


-12.6 d 


8 


2400 s 


1.78 


HIP69366 


1.91 


2011-08-31 


55804.23 


GNIRS 


7.55 d 


-9.7 d 


8 


960 s 


1.77 


HIP69366 


1.77 


2011-09-03 


55807.24 


GNIRS 


10.55 d 


-6.7 d 


8 


1200 s 


1.89 


HIP69366 


2.05 


2011-09-07 


55811.11 


GNIRS 


14.42 d 


-2.8 d 


8 


960 s 


1.88 


HIP69366 


1.99 


2011-09-10 


55814.21 


GNIRS 


17.52 d 


+0.3 d 


8 


960 s 


1.87 


HIP69366 


1.95 


2011-09-13 


55817.24 


GNIRS 


20.55 d 


+3.3 d 


8 


960 s 


2.28 


HIP69366 


2.62 


2011-09-18 


55822.13 


GNIRS 


25.45 d 


+8.2 d 


8 


700 s 


2.25 


HIP69366 


2.54 


2011-09-22 


55826.22 


GNIRS 


29.53 d 


+12.3 d 


8 


960 s 


2.35 


HIP69366 


2.84 


2011-09-27 


55831.21 


GNIRS 


34.52 d 


+17.3 d 


8 


960 s 


2.54 


HIP72520 


2.35 


Note. — 


The time of 


explosion 


of MJD 55796.687 (Nugent et al. 2011]) and the time of B-band maximum of MJD 55813.9 



over six orders. Because of the choice of the l'.'O slit 
width, the resulting spectral resolution is lower than the 
SpeX spectrum and in the range of ~ 1000 — 1300. The 
GNIRS data were calibrated and reduced using the gnirs 
pipeline within the Gemini IRAfF^I package, specifically 
developed for the reduction of GNIRS data. The steps 
began with non-linearity correction, locating the spec- 
tral orders and flat-fielding. Sky subtractions were per- 
formed for each AB pair closest in time, then the 2D 
spectra were stacked. Spatial distortion correction and 
wavelength calibrations were applied before the ID spec- 
trum was extracted. To perform telluric correction, the 
stellar hydrogen lines were first removed from the tel- 
luric star spectrum. The IRAF task telluric was then 
used to interactively adjust the relative wavelength shift 
and flux scale to divide out telluric features present in 
the science spectrum. A blackbody spectrum is then as- 
sumed for the telluric star for the flux calibration. Some 
high-frequency pattern noise with pixel-wide correlations 
were evident in isolated regions of the spectrum; this is 
a known issue for the GNIRS detector. These were re- 
moved with a low-pass filter. There were adequate over- 
laps between orders to join the spectra; however, some 
overlaps occurred in the uncertain regions of strong tel- 
luric absorptions be tween the JHK bands. We used the 
JHK photometry of Mathcson et al. (2012) as a guide to 
scale the relative flux before joining the spectral orders. 

3. NIR SPECTRA OF SN 2011fe 

Ten NIR spectra of SN 201 lfe are presented in Figured] 
H The early NIR spect ra are dominated by electron 
scattering with the continua resembling the Rayleigh- 
Jeans tail of a blackbody. The spectral slope remains 
quite constant until around maximum light. As the 
ejecta expands and electron scattering opacity decreases, 
the continuum drops, and a well-defined photosphere no 
longer exists. Just past maximum light and at regions 
of heavy line blanketing by iron-peak elements, the lines 
are formed at larger e ffective radii produc ing features 
seemingly in emission (IWheeler et alJ fl998l) . The com- 
bination of these lines in "emission," and the decreased 
continuum, forms the most prominent feature in the NIR 

15 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

16 The spectra can be downloaded at http://csp2.lco.cl/hsiao/ 



between 1.5 — 1.9 /im. 

Alternatively, the origin of the increased flux at 1.5 — 
1.9 /im was sh own to be caused by increased fluorescence 
in the region (|Gall et al.ll2012D . Line blanketing of iron- 
group elements is very efficient in blocking flux in the 
ultraviolet (UV) part of the spectrum. In contrast the 
opacity in the optical an d NIR is significantly lower (e.g., 
iPinto &: Eastman! r2000). Thus photons absorbed in the 
UV and redistributed to the NIR by line f l uores cence 
may escape from the ejecta (|Kromer fe Siml [20091 ) . At 
wavelengths where line fluorescence is highly effective, 
this can lead to an increased flux. 

We fit three early NIR spectr a with the automate d 
spectrum synthesis c ode SYNAPPS ( Thomas et al.ll20lTbl) , 
derived from SYN0W ([Branch et al.H2005f ). The results are 



shown in Figure^! SYNAPPS uses a highly parameterized, 
but fast spectrum synthesis technique, useful for identi- 
fying the ions that form the observed features. We in- 
cluded in the fits ions which are commonly identified or 
expected to be present in normal SN la optical and NIR 
spectra: C I, C II, O I, Mg II, Si II, Si III, S II, Ca II, 
Fe II, Fe III, Co II and Ni II. Such an expansive list of 
ions was included, such that SYNAPPS can objectively de- 
termine the absence or presence of each ion species. The 
most prominent features at this epoch are attributed to 
intermediate-mass species, such as O I, Mg II and Si III. 
The strong Mg II A1.0927 /im absorption is formed at a 
relatively isolated region near 1.05 /im. Signatures of un- 
burned C I were detected and will be discussed in detail 
in Section 14.11 Redward of the strong Mg II line near 
1.05 /im, four absorption features from 1.1 /im to 1.4 /im 
were identified as O I/Si III, C I/Ca II, Si III, and Si III, 
respectively. 

In terms of iron-peak elements, SYNAPPS interprets the 
strong feature near 1.6 /xm as Fe HI. This feature was 
identified a s Mg H/Si II/Co II b y iMarion et all (l2009fl 
and Si II bv IWheeler et~all (fl998h and|Gal]_etJI[(j2012). 
Two absorption minima near 1.6 /im were resolved in 
SN 2011fe. Both minima are attributed to Fe III with 
weaker influences from Mg II, Si II and Si III in the 
redder of the two absorptions. The Fe III interpretation 
also yielded good fits for several other features in and 
around the H band (Figure [2]) . Ni II appears to be weak 
in all three early spectra, while Co II increases in strength 
with time. The increasing influence of Co II is attributed 
to both decreasing opacity and increasing abundance of 
56 Co from the radioactive decay of 56 Ni. 
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: expl = 25.45 d 
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;expl 
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= 12.3 d 



2011-09-27 

= 34.52 d 
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Fig. 1. — NIR spectra of SN 2011fe. The wavelength is as observed. The logarithmic scale of flux is plotted. For clarity, the spectra 
are scaled such that they have constant spacing at 1 (im. The UT date of observation, phase relative to explosion and phase relative to 
B-band maximum are labeled for each spectrum. Grey vertical bands mark the regions of the strongest telluric absorptions. 
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The presence or absence of a certain ionization stage 
of a species is regulated by the combination of temper- 
ature and the strength of the lines in question. While 
C II was p ositively identified in the optical spectra of 
SN 201 lfe (jParrent et al.ll2012D , the C II lines are weak 
in the NIR. Instead, C I is identified in the NIR, with 
the strongest influence near maximum light rather than 
at earlier phases (Section 14. ip . While Si II is the hall- 
mark feature of a SN la in the optical, Si III was found 
to have more identifiable features in the NIR spectra of 
SN 201 lfe. At early phases, Fe III has stronger influ- 
ences than Fe II in the NIR, similar to the behavior in 
the optical. 

Recently, iGall et all (|2012f) showed the close resem- 
blance between the NIR spectra of SN 2005cf and the 
synt hetic spectra of the hydrodynamical explosio n model 
W7 (|Nomoto et all H981 llwamoto et all 11999( 1. They 
found that in the NIR, the emerging flux is almost en- 
tirely due to fluorescent emission with a few features 
formed by strong P-Cygni absorptions. Here, we make 
the comparison with the NIR spectra of SN 201 lfe, using 
the same synthetic spectra without modifications (Fig- 
ure [3|). While SYNAPPS is useful for identifying the line 
forming ion species, it does not provide detailed post- 
processing of stellar explosion models. The W7 model 
spectra in Figure[3]were calcula ted using the Monte C arlo 
radiative transfer code ARTIS (jKromer fe Sim! 12009ft . It 
provides a fully self-consistent solution to the radiative 
transfer problem and produces time-dependent spectral 
synthesis without any free parameters. As was found by 
IGall et aT] (|2012f ) . the spectral features of the observed 
and model spectra show encouraging similarities. A few 
differences are noted here. Slight discrepancy is observed 
in the time evolution of the NIR brightness. The two 
spectra of SN 201 lfe near a month past explosion or 10 
days after B-band maximum are consistently fainter than 
the model spectra. The strong Mg II A1.0927 /im fea- 
ture near 1.05 /im is present in the early model spectra; 
however, discrepancies in the velocity evolution are evi- 
dent between the observed and model spectra. While the 
Mg II A1.0927 /im velocity of SN 2011fe decreases rapidly 
(Section 14. 2[) . the velocity in the model spectra remains 
high. Also, the observed strong break near 1.5 fim (Sec- 
tion !4.3[) . one of the most prominent spectral features in 
the NIR, is not reproduced by the model until approx- 
imately one month past the explosion. Even then, the 
velocity of the feature does not match the observations. 

4. DIAGNOSTICS OF SN la PHYSICS 

There is currently no full theoretical description of 
the explosion mechanism of SNe la that accounts 
for the observed diversity. In the most established 
model, carbon is ignited near the center of a carbon- 
oxygen i _white_dwarf ; _as it app roaches the Chandrasekhar 
mass (Hoylc fe Fowlerl 11960ft . The initial thermonu- 
clear carbon burning starts in a subsonic deflagration 
(|Nomoto et alj 119761 ). To produce partially-processed, 
intermediate-mass material observed in SNe la, the tran- 
sition to detonation is "delayed," such that the white 
dwarf is able to expand appreciably to allow for partial 
burning. A possible transition to a su personic detona- 
tion then rapidly inc inerates the star (|Khokhlovl 119911 : 
lYamaokaet aH 119921) . While some mechanisms for the 
transition from deflagration in the carbon flame to the 



subsequent detonation in the oxygen flame have been 
prop osed, its nature remain s uncertain and controversia l 
(e.g.. iKhokhlov et all 119971; iNiemever fe Wooslevl Il997t 
Plewa et all 120041; lAspden et all I2008t IWooslev et all 
2009, 2011). The conditions under which the transi- 



tion occurs are still essentially treated as free parameters 
within models. 

Early NIR spectra provide several powerful diagnos- 
tics of the physics of SN la explosions, and can serve 
as discriminators for competing explosion models. Since 
magnesium is a product of explosive carbon burning and 
not oxygen burning, the NIR Mg II lines provide a di- 
rect measure of the inner boundary of carbon burning, 
and could place meaningful constraints on the condi- 
tions under which th e transition to detonation occurs 
(| Wheeler et all 119981 ). The location of the boundary 
is also particularly sensitive to the transition density, 
which is ine xtricably linked to the am ount of 56 Ni pro- 
duced (e.g.. IHoflicheTail I1995L 120021) . The presence or 
absence of unprocessed material is one of the key pre- 
dictions of SN la explosion models. The NIR C I lines 
provide an independent confirmation of the findings from 
the intensively-studied optical C II lines. Past maximum 
light, the continuum opacity from electron scattering de- 
creases dramatically. Coupled with heavy line blanketing 
from fully-processed, iron-peak material, the NIR spec- 
trum probes a large range of de pths and the tempera ture 
of the line blanketing material jWheclcr et al. 1993). In 
this section, we present several techniques to quantify 
the characteristics of NIR spectra and discuss the use 
of these measurements as various diagnostics of SN la 
physics. 

4.1. Unburned Carbon 

Since oxygen is also produced from carbon burning, 
carbon provides the most direct probe of the primor- 
dial material from the progenitor carbon-oxygen white 
dwarf. The quantity, distribution and incidence of un- 
burned carbon in SNe la provide important constraints 
for explosion models. Turbulent deflagration models pre- 
dict that a l arge amount of unprocessed carbon s hould 
be left over (|Gamezo et alj|2003t iRopke et al.ll2007D . In 
general, a transition t o detonation would result in com- 
plete carbon burning (Hofli ch et all [2002; M arion et ail 
l2006t [Kasen et al. 2009); although, exceptions have been 
found. For the lowest transition density in a grid of de- 
layed detonation models, both the model and the ob- 
servations of SN 19 99by showed strong C I signatures 
(iHoflich et al.ll2002l). but SN 1 999by is significantly sub- 
luminous dToth fe S zabo 2000) a nd n ot representative of 
a normal SN la. iGamezo et alj (|2004ft showed that pock- 
ets of cool and dense unburned material could be pulled 
down toward the center and escape the detonation wave. 

Unburned carbon in SNe la has been the subject of 
several intensive studies in recent years. Examining the 
C I lines in NIR spectra of three normal S Ne la without 
the ai d of spectrum synthesis techniques, iMarion et al.l 
(2006) conclud ed that the ab u ndanc e of unprocessed ma- 
terial is low. [Thomas et al.l (|2007D presented the most 
convincing detection of the C II AO. 6580 /im feature 
in the early optical spectra of SN 2006D. In subse- 
quent studies with larger optical spectroscopic samples, 
20 — 30% o f the early spectra we r e found to show C II 
signatures (|Thomas et al.l l2011al : iParrent et al.l 120111 ; 
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Fig. 2. — SYNAPPS fits for three early NIR spectra of SN 2011fe. The observed spectra and the SYNAPPS fits are plotted in black solid and 
grey dotted curves, respectively. The colored curves show individual contributions from the most prominent ions. The wavelength axis is 
plotted in the rest frame of the host galaxy. The phase relative to B-band maximum is labeled for each spectrum. 



IFolatelli et aHl20ll [Silverman fc Filippenkoll2012l) . Sev- 
eral studies have noted that this fraction represents the 
lower limit, as noise, velocity blueshift, and the phase 
at which the SNe la were followed up could affect the 
C II detection (e.g., iBranch et all l2007t IFolatelli etall 
2012). The velocities of the detected C II lines are gen- 
erally low. There have also been hints that SNe la with 
detected C II have preferentially bluer colors and nar- 
rower light curves {Th omas et all 1201 lal IFolatelli et all 
20121). but this finding was not h eld up in the exam- 
inations of other data sets (e.g., iBlondin et al.l 120121 : 
iParrent et al.ll20lTI ). Nonetheless, there appears to be a 
consensus that t he mass fraction of the photospheric car- 
bon is low (e.g.. IBranch et al.1120031: iMarion et al.l l2006t 
iThomas et al.l 120071 : iTanaka et al.ll200l ~ 

The automated spectrum synthesis code SYNAPPS has 
been employed successfully to iden tify the C II A0.6580 
/im lines in early op tical spectra (| Thomas et al.1 120071 
1201 lal IParrent et alJl20ll . Here we concentrate on the 
study of the strongest C I line in the NIR at 1.0693 /im 
using SYNAPPS. The weaker C I lines are included in the 
fit, but are heavily blended and more difficult to detect. 
For example, the C I AA0. 9087, 1.1756 /im lines are both 
blended with O I, Si III, and Ca II lines. C I A1.0693 /im, 
at the largest expected range of velocities, lies in between 



the strong Mg II A1.0927 /mi line, usually blueshifted to 
~ 1.05 /tm, and the weaker Mg II A1.0092 /im line, usu- 
ally blueshifted to ~ 0.97 /im. The fits were done for 
three early spectra over the entire NIR wavelength range 
(Figure [2j . The results zoomed in on the region con- 
taining the line profiles of C I A1.0693 /im are shown 
in Figure 2] The best-fit synthesized spectra show very 
close matches to the observed spectra. In the first spec- 
trum, at approximately 2 weeks before maximum, the C I 
signature appears weak, and was not required to obtain 
a decent fit in the region. Starting from a week before 
maximum, the story changes. The profile begins to take 
on a flattened shape (Figure 13}. The presence of C I is 
now required to suppress the emission component of the 
Mg II A1.0092 /im P-Cygni profile, which is evident at 2 
weeks before maximum. Near maximum light, the pro- 
file remains flat, and the influence of C I is even stronger. 
Our spectrum has a high enough signal-to-noise ratio to 
resolve the "notch" from C I A1.0693 /tm near 1.03 /tm. 
The strengthening of the C I signature with time may be 
a result of the increasing ionization fraction of C I ver- 
sus C II, as the ejecta expands and cools. The unburned 
carbon in these early NIR spectra was detected at veloci- 
ties consistent with the photospheric velocity of partially 
and completely synthesized material. The blueshift of 
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Fig. 3. — Comparison between the synthetic W7 model spectra 
and the NIR spectra of SN 2011fe. The model spectra were cal- 
culated using the radiative transfer code ARTIS for the W7 model. 
The observed and the model spectra are plotted in black and red, 
respectively. The wavelength axis is plotted in the rest frame of the 
host galaxy of SN 2011fe. The absolute flux scales of the observed 
and model time series were matched at the first comparison epoch, 
then the same offset was applied to each subsequent pair of spectra. 
The phase relative to explosion is labeled for each spectrum. 

the C I line was observed at 12,900 kms^ 1 at 2 weeks 
before maximum and leveled off at ^12,000 kms -1 from 
a week before maximum light. Note that in past analyses, 
the flattened feature has been attributed to an ex tended 
wing of Mg II Al. 0927 /mi (e.g. JMarion et al.|[2009l ). This 
is produced by a detonation front which leads to parti al 
burning up to very high velocities ; I loliich et al.l l2002). 

Our fits were performed without assuming any initial 
conditions based on the results from the optical, such as 
strength and velocity. In Figure the C I A1.0693 /im 
feature from the NIR and C II AO. 6580 /im feature from 
the optical of SN 2011fe are compared in velocity space. 
Optical spectra taken at phases closest to those of the 
NIR spectra are shown. The C II features were posi- 
tively identified in the early optical spectr a of SN 20 life 
and s howed decreasing strength with time (jParrent et al.l 
2012). The velocities of the C I and C II features match 
exceptionally well, indicating that these signatures are 
from the same line-forming material. The interpreta- 
tion from SYNAPPS that C I A1.0693 /im is the source of 
the flattened profile on the emission wing of the Mg II 
A1.0092 /im P-Cygni profile appears to be correct. 

It is worth noting here that the flattened profile near 
1.03 /im is a common feature for brighter and slower- 
declining SNe la near maximum light (Figure |5J). For 
fainter, faster-declining SNe la, the C I signatures ap- 
pear to be stronger and easily detectable at maximum 
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Fig. 4.— SYNAPPS fits of three early NIR spectra of SN 2011fe in 
the region of the C I A1.0693 fim line. The velocity axis is plotted 
with respect to the rest wavelength of C I A1.0693 fim. The spectra 
are plotted as solid black curves. The best-fit synthesized spectra 
are plotted as follows: with all ions, with only C I, and with all ions 
except C I. These are plotted as red dotted, green dashed, and blue 
dash-dotted curves, respectively. The vertical dotted lines mark 
the location of the best-fit minimum photosphcric velocity. The 
phases relative to ZJ-band maximum light are noted. 

light. SN 1999by is the extreme of this example with 
conspicuous presence of the C I line. This effect is rem- 
iniscent o f the familiar seque nce of the Si II AO. 5972 /im 
strength ([Nugent et al.lll995h . It is possible to interpret 
from these observations that the majority of SNe la, at 
both ends of the brightness extreme, harbor some un- 
processed carbon deep in the ejecta. Such a conclusion 
is discrepant from the 20 — 30% fo und by recent optical 
studies (e.g.. iThomas et al.ll2011at IFolatelli et"aLl l2012t 
iSilverman fc FilippenkdT2^12T )7 although, there are indi- 
cations that the fraction of detection is higher with ear- 
lier spectra (e.g.. IFolatelli et aT]|2012D . Since the finding 
of ubiquitous unprocessed carbon would have profound 
implications for our understanding of SN la explosions, 
a larger NIR sample and further studies are warranted 
for such a claim. 

Nonetheless, we have shown here that the C I A1.0693 
/xm feature in the NIR is potentially a better tracer of 
unprocessed material than C II AO. 6580 /im in the op- 
tical. While the signature of C II AO. 6580 /im is the 
strongest between one and two weeks before maximum 
and fades with time, the influence of C I A1.0693 /im ap- 
pears to be the strongest at maximum light in SN 2011fe. 
A change in the ionization condition, as the temperature 
cools, may have fortuitously brought about a delay in 
the onset of the C I feature. If this behavior can be veri- 
fied in other SNe la, the problem of the incompleteness of 
very early optical spectroscopic samples could be circum- 
vented by using a sample of NIR spectra near maximum 
light. Well situated in a relatively isolated region be- 
tween the two Mg II lines, the NIR C I line also has the 
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Fig. 5. — Comparison of the NIR C I A1.0693 fim and the optical 
C II AO. 6580 (im lines of SN 2011fe in velocity space. The velocity 
axis is plotted with respect to the rest wavelength of the C I A1.0693 
fim and C II AO. 6580 fim lines. The NIR spectra are plotted as 
black curves. The optical spectra are taken from IParrent et alj 
(2012) and plotted as red curves. The C I profiles from the best- 
fit synthesized spectra are plotted as black dashed curves. The 
vertical dotted lines mark the location of the best-fit minimum 
photospheric velocity. The phases relative to B-band maximum 
light for the NIR and optical spectra are labeled in black and red, 
respectively. 

advantage that it is unblended with other strong absorp- 
tions for a large range of velocities; while C II AO. 6580 
fim may be blueshifted into Si II AO. 6355 fim at a velocity 
higher than 15,000 kms" 1 (|Branch et al.ll2007ft . 

4.2. Inner Edge of Magnesium 

IWheeler et al.l (|1998D identified the strong and rela- 
tively isolated absorption feature at 1.05 fim as Mg II 
A1.0927 fim. The line is expected to be observed with de- 
creasing velocity. The velocity is then expected to cease 
changing when the photosphere has receded below the in- 
ner edge of the Mg II distribution. Magnesium is a prod- 
uct of explosive carbon burning, but not oxygen burning, 
as the required higher temperatures for oxygen burning 
generate products further along the alpha chain, such as, 
silicon and sulfur. This makes the Mg II line a sensitive 
probe to the location of the inner edge of carbon burn- 
ing in velocity space (IWheeler et al.l 1998; M arion et al.1 
120011: iHoflich et ail 120021 : iMarion et al.ll2009f ). Here, we 
focus on measuring the velocity using the location of the 
absorption minimum of the Mg II A1.0927 fim line pro- 
file. With the NIR spectra of SN 201 lfe and previous 
data from the literature, we explore the range of Mg II 
velocity and the accuracy limit of the velocity measure- 
ments. 

The canonical Gaussian line profile is inadequate in 
fitting the Mg II P-Cygni profile, which contains influ- 
ences from the P-Cygni wings of C I, O I and Si II. Using 
the Gaussian function to fit supernova line profiles often 
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Fig. 6.— Comparison of the NIR C I A1.0693 ,m profile for 
five SNe la near maximum light. The velocity axis is plotted with 
respect to the rest wavelength of C I A1.0693 fim. The C I profile 
from the best- fit synthesized spectrum of SN 201 lfe is plotted as 
a black dashed curve. The time of B-band maximum and B-band 
light-curve decli ne rate Amisf B) of the literature spectra were 
taken from Hickcn et al. (2009). The ph ase relative to B-band 
maximum and Ami5(B) are noted for each spectrum. 



results in poor model fits which are unaccounted for by 
the errors. To properly estimate the uncertainties in the 
velocity measurements, we opted to use a principal com- 
ponent model which was built from data and empirically 
describes the variation in the line profiles. The formula- 
tion and the construction of the model are described in 
the appendix. 

The first four principal components and the model 
spectra, which account for over 90% of the variations, 
are plotted in Figure [7J The first two principal compo- 
nents, £o and £i, show strong correlations between the 
strengths of the Mg II AA1.0092, 1.0927 fim lines and 
the C I A1.0693 line. These line strengths are also 
correlated with those of O I and Si III lines further to 
the red. The next principal component, £2, describes ap- 
proximately the varying tilts in the input spectra. Using 
these principal components, the model spectrum of the 
Mg II feature was then built as M(p,n), which varies 
as a function of the number of principal components n 
and the projections p onto the corresponding principal 
components. 

To fit the observed spectrum, denoted by O(A) and 
with uncertainty AO(A) and m wavelength elements, we 
allowed extra freedom in the model spectra M(n,p) in 
the normalization C and wavelength shift A s hift , and min- 
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Fig. 7. — Principal component model of the region centered on 
the Mg II A1.0927 /jm line. The first four principal components 
are shown. For the model spectra, the effects of 1 <r variation in 
each projection pj are shown. The eigenvalue for each principal 
component is noted on the top and is expressed as a fractional 
value, such that the sum of the N eigenvalues totals to 1. It is 
a measure of the amount of variation each principal component 
describes. 
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The best fit model spectrum was det ermined using the 
IDL package mpf it (Markwardt 2009), a nonlinear least- 
squares fitting program ported from the FORTRAN 
package MINPACK-1. The number of principal compo- 
nents to include, n, was fixed in the fit and chosen, such 
that, the reduced x 2 of the resulting best fit is closest to 
1, typically n = 6 — 8. 

Once the best-fit model spectrum was found, we at- 
tempted to remove the steep "co ntinuum" with meth - 
ods similar to that employed by iMarion et al.l (|2009f ) . 
A straight line connecting the model spectrum at fixed 
wavelengths to the blue and red side of best-fit A s hift 
was taken as the continuum and removed from the 
model spectrum. The profile minimum of the continuum- 
removed model spectrum was then used to determine the 
line velocity of Mg II. The fit uncertainty of A s hift was 
used to derive the measurement error for the velocity. 

The Mg II velocities from several literature early-phase 
time-series NIR spectra were determined and plotted 
wit h the velocities of SN 2011fe in FigureEl As predicted 
bv IWheeler et al.l (119981 ). SN 201 lfe shows a rapid de- 
crease in the Mg II line velocity followed by an extended 
epoch of constant velocity, beginning at ~ 10 days be- 
fore maximum and lasting until the feature disappears 
at ~ 10 days past maximum. The velocity during this 
phase is remarkably constant for SN 201 lfe, showing a 
dispersion of only 130 kms -1 , on the order of the veloc- 
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The time evolution of the Mg II A1.0927 fim velocity. 



ity resolution of these medium-resolution spectra. The 
emission portion of the C I A1.0693 /im P-Cygni profile is 
blended with the absorption feature of Mg II A1.0927 [im 
(Section |4Tj. In the case of SN 201 lfe, the time- varying 
characteristic of the C I feature had very little influence 
on the Mg II velocity during this epoch. For the rest of 
the SNe la in our sample, they appear to be caught after 
the Mg II line has entered the epoch of constant velocity, 
with the exception of the subluminous SN 1991bg-like 
SN 1999by. For SN 1999by, the Mg II velocity contin- 
ues to decrease well past maximum light and indicates 
a distribution of Mg II down to ~ 8,000 kms -1 . Note 
that SN 2002dj is a significant outlier in Figure [5] with 
exceptionally high Mg II velocity which persisted for a 
wide phase range before maximum light. SN 2002dj is 
also identified as a high velocity g radient object in th e 
IBenetti et all (j2005[ ) cl assification (IPignata et al.H 2008D, 
similar to SN 2002bo (jBenetti et al.l I2004Q . The NIR 
spectrum of SN 2002bo also shows high Mg II velocity, 
although not as extreme as that of SN 2002dj. These 
results may point to a correspondence between high Si II 
velocity gradients and high Mg II velocities. 

Here, we test whether the location of the inner edge of 
explosive carbon burning has a significant impac t on the 
56 Ni productions of SNe la. Using Ami 5 (B) (jPhillips! 
[199! as a proxy for the amount of 56 Ni produced, the 
result is plotted in Figure [9l The Mg II velocity for each 
SN is the mean velocity, weighted by the measurement 
errors, in its constant epoch. The uncertainty is the in- 
trinsic dispersion of the velocity around the mean during 
the constant epoch. For those SNe without adequate 
time-series data, the uncertainty is assumed to be the 
mean of the dispersions of all SNe la with time series, 
excluding SN 1999by. 

There is a large spread in the locations of the inner 
edge of the Mg II distribution in SNe la. For a range of 
normal SNe la 1.0 <Arni 5 (B)< 1.2, the spread of Mg II 
vel ocities is as h i gh as 5,000 kms -1 , as was observed 
by IMarion et all (|2003l 12009ft . There appears to be no 
correlation between the Mg II velocity and Ami5(B). 
This result is surprising, since the transition density is 
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Fig. 9. — The Mg II A1.0927 fim velocity versus optical light- 
curve decline rate Ami5(_B). SNe la with three or more NIR spec- 
tra taken during the constant velocity epoch are plotted in filled 
circles, while the rest are plotted in open circles. For the peculiar 
SN 1999by, a downward arrow indicates that the constant veloc- 
ity epoch has yet to be reached in the defined time interval (see 
Figure [HJ. 



expec ted to have a strong influence on the yield of 56 Ni 
( [ e J g ^Hoflich"et~aT] fl995l) . The model of IHoflich etafl 
(2002) for SN 1999by corresponds to the lowest transi- 
tion density in the model grid. SN 1999by also happens 
to be a subluminous event and has very low Mg II ve- 
locities. Yet SN 2005am, with the next fastest decline 
rate in our sample, has one of the highest Mg II veloc- 
ities measured. If the Mg II velocity is indeed a strong 
function of the transition density, its effect on the 56 Ni 
production and the brightness of SNe la appears to be 
secondary. Since its correlation to Ami5(-B) is weak, it 
provides independent information which may affect the 
brightness of SNe la. We explore the possibility of using 
the Mg II velocity measurement to further improve the 
use of SNe la as distance indicators in Section [6] 

4.3. H-band Iron-peak Feature 

Soon past -B-band maximum light, the NIR spec- 
tra take on a characteristic shape. Deep depression 
in regions between 1.2 — 1.5 /tm is coupled with a 
dramatic feature between 1.5 — 1.9 /an, seemingly in 
emission (Figure [J). The feature was first noted by 
IKirshner et all (|l973h . We now understand the deep de- 
gression as the result o f a lack of line-blanketing opacity 
Spyromilio et al.ll994D . The early NIR spectra of SNe la 
are dominated by electron scattering. As the ejecta ex- 
pand and cool, the electron scattering opacity decreases 
and, in turn, lowers the continuum flux level and exposes 
the core. In contrast, the strong line blanketing opacity 
from iron-peak elements causes the prominent feature at 
1.5 — 1.9 /mi to fo rm at a larger effec tive radius and 
appear m emission (|Wheeler et al.lll998D . The strongly- 
variable opacity provides probes of very different depths 
at the same epoch. An alternate origin of the increased 
flux at 1.5 — 1.9 fim was proposed to be i ncreased fluo- 
rescence in the region bv lGall et al.l (|2012[) . 
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Fig. 10. — Time evolution of the ratio across the break at 1.5 
(im R. The inset in the top right corner illustrates the definition 
of the break ratio with the spectrum SN 2011fe at 12.3 days past 
maximum. 



In the time-series data of SN 201 lfe, we caught the 
budding of the prominent iron-peak feature in the spec- 
trum at ~ 3 days past B-band maximum (Figure [1]). It 
then emerges rapidly to be the dominant feature in the 
NIR. This feature coincides in wavelength with the H 
band. The rapid time evolution and the shear size of 
the feature inevitably cause large uncertainties in the H- 
band k-corrections. Understanding the evolution and the 
diversity of the feature is crucial to control the system- 
atic errors originating from k-correction estimates. The 
relative flux levels in regions of contrasting opacity also 
provide a temperature diagnostic for the line-blanketing 
material. Here, we explore the flux ratio across the break 
at 1.5 /im, which separates the different opacity regimes. 

We define the ii-band break ratio R as the ratio of 
the maximum flux just redward of the 1.5 /jm break to 
the minimum flux just blueward of the 1.5 /im break 
(Figure [TU)) . To determine the flux ratios for the sample 
of NIR spectra, the spectra were first smoothed with a 
Gaussian filter with a fixed width of AA/A = 0.0025 in 
logarithmic wavelength space. The minima and max- 
ima were found in fixed wavelength regions from the 
smoothed spectra. Since the velocity of the feature and 
the associated uncertainty do not enter in this analysis, 
we found that the Gaussian smoothing technique to be 
simpler and more robust than the PCA models described 
in Section 14.21 The measurement errors were estimated 
from the dispersions in the neighboring pixels, within 
±100A, of the minima and maxima. The signal-to-noise 
ratio of each spectrum around the H-band break region 
is therefore reflected in the error. The resulting if-band 
break ratios for SN 2011fe and the sample of literature 
spectra are plotted in Figure [TU] The ff-band features 
of SNe la appear to emerge at very similar phases, ~ 3 
days past -B-band maximum light. With the exception 
of the peculiar SN 1999by, the flux ratios appear to peak 
uniformly at ~ 12 days past B-band maximum light; 
although, more time-series data are required to charac- 
terize the rise. Then, a remarkably uniform linear decline 
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Fig. 11. — The ratio across the break at 1.5 fim at 12 days past 
maximum, R12, versus optical light-curve decline rate Ami^(B). 
SNe la with two or more time-series observations in the post-peak 
decline are plotted with filled circles, while the rest are plotted with 
open circles. 



is observed among normal SNe la. Excluding the pecu- 
liar SN 1999by, we determined an error- weighted mean 
linear decline of —0.146 day -1 with a standard deviation 
of 0.036 day" 1 . 

From our limited sample, the time evolution of the H- 
band break ratio R appears to be uniform with a range 
of values at peak. We defined R12 as the -ff-band break 
ratio at 12 days past -B-band maximum. For each SN la 
with time-series data after the ratio peak, a straight line 
was fit to the linear post-peak decline. Then R12 was 
determined by extrapolating on the best-fit line to 12 
days past B-band maximum. The uncertainty in R12 
was derived from the error in the slope of the linear fit. 
For SNe la without adequate time-series data to deter- 
mine the linear decline of the if-band break ratio, the 
error-weighted mean decline rate of —0.146 day^ 1 and 
an average uncertainty were adopted. The resulting R12 
measurements are plotted against Amis (2?) m Figure [TO 

The peak iJ-band Break ratio shows strong correla- 
tion with Ami5(£?) (Figure [TT]): although, more time- 
series data are needed to confirm this finding. If the 
effect from intermediate-mass material at 1.2 — 1.5 /xm 
is secondary, the correlation describes a temperature se- 
quence. The temperature of the iron-peak material ap- 
pears to be closely tied to the production of 56 Ni, more 
so than the transition density (Section I4.2[) . In the view 
point of contrasting opacity in the region, the observed 
effect can be explained as follows. At ~ 12 days past 
B-band maximum, the electron scattering opacity is low 
enough that the iron-peak feature is fully exposed. The 
temperature of the iron-peak material now dictates the 
size of this feature. As the ejecta expands, the line- 
blanketing opacity drops with the decreasing tempera- 
ture of the iron-peak material. The iron-peak feature at 
1.5 — 1.9 fim then forms at decreasing effective radii and 
appears weaker. The uniformity of the time evolution 
indicates that the temperature decline of the iron-peak 
material is consistent in normal SNe la, regardless of the 
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Fig. 12. — Comparison betwe en the NIR spe ctra of SN 2011fc 
and the revised template of Hsiao ct al. (2007) warped to match 
observed colors from izYJHK photometry. At the bottom, the 
locations of the observed NIR filters at z = and z = 0.08 are 
plotted in the rest frame of the supernova. 

varying amounts of 56 Ni produced. The predictability 
of the strength and the evolution of the iron-peak fea- 
ture means that we can indeed improve upon the c urrent 
meth od of NIR k-correction (jHsiao et al.l 120071 iHsiaol 
2009), by building, for example, spectral templates as 
a function of Amis(i3). 

5. NIR K-CORRECTIONS 

K-corrections account for the effect of cosmological ex- 
pansion. SNe la at different redshifts would have their 
spectral energy distribution (SED) sampled at different 
wavelength regions by the same observed filter. A k- 
correction converts an observed magnitude to that which 
would be observed in the rest frame of the same or an- 
other filter, allowing for the co mparison of the brightn ess 
of SNe la at various redshifts (|Oke fc Sandaed figesFI. 
The calculation of k-corrections requires knowledge of the 
SED. While the effect of an assumed spe ctral template 
has been extensively stu died in the optical (|Nugent et al.l 
120021 iHsiao et al1l2007j ). the effect is not well understood 
in t he NIR due t o the relatively small spectroscopic sam- 
ple (|Hsiaoll2009| ). 

The current method of NI R k-corrections em ploys the 
revised spectral templates of iHsiao et al.l ()2007f ) , warped 

17 See Aldcrine ct al. (2002) for an alternative experimental ap- 
proach which eliminates the need for k-corrections. 
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Fig. 13.— NIR YJHK k-c orrection difference between SN 2011fc 
and the revised template of Hsiao et al. (2007). Two redshifts are 
shown: z = 0.03 (filled circles and solid lines) and z = 0.08 (open 
squares and dotted lines) . The horizontal lines mark the estimated 
offsets in the peak magnitude caused by the k-correction difference 
for z = 0.03 (solid lines) and z = 0.08 (dotted lines). 



to m atch the observed colors of the SN la ([Burns et al.l 
120111 ). Here, we explore the effect of assuming the spec- 
tral template for the case of a SN 2011fe-like object at 
various redshifts. Before the k-corrections were calcu- 
lated, the spectral template at each phase was warped 
to the NIR colors provided by WFCAM and WHIRC 
izYJHK photometry of SN 2011fe (Im et al. in prep; 
iMatheson et~aH 120121 ). The NIR spectra of SN 2011fe 
and the color-corrected spectral templates are plotted 
together for comparison in Figure 1121 

The mismatch of the Mg II feature near 1.05 jxm is 
evident in the first three spectra. The discrepancy is 
not surprising, since these are the earliest NIR spec- 
tra published to date. There also appears to be some 
mismatch in the strength of the prominent iJ-band fea- 
ture at later phases. This illustrates the limitation of 
using a single template for all SNe la. As shown in Sec- 
tion 14.31 the strength of this iron-peak feature correlates 
with Arni5(B). Armed with this well-defined correla- 
tion, one can build a sequence of templates as a function 
of Ami§(B) to minimize the k-correction uncertainty as- 
sociated with this feature. 

The mismatches between the spectra and template are 
reflected in the k-correction differences plotted in Fig- 
ure [T31 The k-corrcction differences were calculated at 
two redshifts, 0.03 and 0.08, which coincide with the red- 
shift limits of current surveys of nearby SNe la in the 
Hubble flow. The larger differences observed at early 
phases in the Y band and at late phases in the H band 
originate from the mismatches in the Mg II feature at 
early phases and the iron-peak feature at late phases. At 
z = 0.08 where the rest and observed filters are signif- 
icantly misaligned, the mean absolute sizes of the dif- 
ferences over the phase range are 0.04, 0.03, 0.04, 0.06 
mag in Y, J, H, K, respectively , similar to the num- 
bers reported bv lHsiao et al.l ([2007D for the optical filters. 
The offsets in the peak magnitude caused by assuming 



the spectral template can be approximated by using the 
average of the k-correction differences weighted by the 
light-curve template errors. For a SN 2011fe-like object, 
the impact of assuming a template SED on the determi- 
nation of peak magnitude is small, on the level of 0.01 
mag for z = 0.03 and 0.02-0.04 for z = 0.08 (Figure [TBI) . 

Nonetheless, whether the differences at various phases 
correlate with each other offers us important clues as to 
how to proceed to minimize the impact of uncertainties in 
k-corrections. If the k-correction differences from various 
phases are completely uncorrelated, more observations 
on the light curve could drive down the uncertainty in 
the peak magnitude determination. If the differences are 
significantly correlated, there is room for improvement 
in the time evolution of the spectral templates. The first 
attempt to quantify these c orrelations in the optical was 
made bv lHsiao et al.l ([2007D . This study requires a large 
sample of time-series spectra, which does not yet exist in 
the NIR. In Section [51 we describe the ongoing efforts to 
build such a data set. 

The availability of color constraints from neighboring 
filters is also crucial for obtaining good k-correction es- 
timates. For example, an accurate k-correction to the J 
band requires Y — J and J — H colors to warp the spec- 
tral template to the correct observed colors around the 
wavelength region in question. Missing color informa- 
tion inflates the k -correction uncertainty to ~ 0.1 mag 
([Hsiao et al.l 12007}) . The higher k-correction difference 
seen in the K band in Figure [T31 is most likely due to the 
lack of constraining observation in a redder band. 

6. CONCLUSIONS 

We have presented ten medium-resolution, high signal- 
to-noise ratio NIR spectra of SN 2011fe. This data set 
constitutes the earliest time-series NIR spectroscopy of a 
SN la, with the first spectrum obtained at only 2.58 days 
past the explosion. We take advantage of these densely- 
sampled observations to gain insights into the time evo- 
lution and diversity of SN la NIR spectral features. The 
main findings are outlined below: 

1. With the aid of the automated spectrum synthe- 
sis code SYNAPPS, C I A1.0693 ^m is detected 
with increasing strength up to maximum light in 
SN 2011fe. The delay in its onset demonstrates 
the potential of the NIR C I line to be a superior 
tracer of unprocessed material to the more com- 
monly used optical C II AO. 6580 fim line. 

2. The Mg II velocity of SN 2011fe, measured from the 
absorption minimum of the Mg II A1.0927 /jm line 
profile, decreases rapidly at very early phases, then 
flattens to a remarkably constant evolution from 10 
days before to 10 days past B-band maximum, with 
a dispersion of only 130 kms -1 . 

3. Taking advantage of the constant velocity over a 
wide time window, the Mg II A1.0927 ^m line, even 
without time-series observations, can be used to 
locate the inner edge of carbon burning. The Mg II 
profiles of SNe la exhibit a wide spread of velocities, 
as high as 5,000 kms -1 , and show no correlation 
with Ami5(B). 

4. The prominent feature at 1.5 — 1.9 fim shows uni- 
form phase evolution in SNe la. The flux ratio 



NIR spectroscopy of SN 2011fe 



13 



across the break at 1.5 /zm increases from ~ 3 
days past maximum, when the break appears, and 
reaches its peak at ~ 12 days past maximum. As 
a proxy for the temperature of the line-blanketing 
material, the strong correlation between the peak 
flux ratio and Ami5(£?) indicates a close tie be- 
tween temperature of the iron-peak elements and 
the production of 56 Ni. 

5. Despite the presence of the prominent and rapidly- 
evolving feature at 1.5 — 1.9 fim near maximum 
light, the predictability of its strength and the evo- 
lution means that NIR k-corrections can be im- 
proved by producing spectral templates as a func- 
tion of Amis(-B). 

The result that the Mg II velocity does not correlate 
with A?ni5(i?) is an intriguing one, since the transi- 
tion density is expecte d to have a stron g influence on 
the yield of 56 Ni (e.g., IHoflich et all 119951) . Is it possi- 
ble that the Mg II velocity describes variations in SN la 
luminosity on a secondary level? Recently, SN la lumi- 
nosity has been found to have a dependence on the host- 
galaxy stellar mass, in addition to the well-k nown depen- 
dence on light-curve shape and colors (e.g.. iKellv et al.l 



2010; Sulliva n"etalll2OT0l: lEarnpeitl et al 



120101) . It has 

been postulated that the effect may be attributed to 
the i mpact of progenitor metallicity on 56 Ni pr oduction 
(e.g., iTimmes et all [20031: Uackson et all l2010h . With 
the transition densi ty linked to progenitor metallicity 
(| Jackson et all 120101) . the NIR Mg II velocity is poten- 
tially a more direct probe of progenitor metallicity and 
a better calibrator of SN la luminosity than host-galaxy 
stellar mass. 

With the limited size of the current NIR spectroscopic 
sample, the Mg II velocity shows a weak trend with the 
host-galaxy stellar mass, hinting a link of Mg II velocity 
to progenitor properties; although, a larger sample with 
unbiased host-galaxy properties is needed to determine 
if such a correlation indeed exists. To test whether the 
Mg II velocity provides better calibration of SN la lumi- 
nosity than host-galaxy stellar mass, a statistical sample 
of NIR spectra, especially those of SNe la in the Hub- 
ble flow, is required. With improved NIR spectrographs 
available on large telescopes, obtaining NIR spectra of 
SNe la in the Hubble flow is now feasible. For example, 
a pre-maximum NIR spectrum of SN 2012ar, a SN la 
with a recession velocity of 8472 kms -1 , was obtained 
in the low-resolution mode with the Folded-port Infrared 
Echellette (FIRE) on the 6.5-meter Magellan Telescope. 
The spectrum reaches a signal-to-noise ratio of > 50 at 
1 /an, with an on-target integration time of only 10 min- 
utes. With this signal-to-noise ratio, the measurement 
error for the Mg II velocity is still well below the in- 
trinsic velocity dispersion of SN 20 life in the constant 



velocity epoch (Section [ 

Throughout the paper, we have emphasized the need 
for an improved data set of NIR spectra to study this rel- 
atively unexplored area and to confirm the various find- 
ings from this paper. Starting in late 2011, we began a 
five-year program to build a large sample of high qual- 
ity SN la NIR spectra using NIR spectrographs on large 
telescopes, such as FIRE, GNIRS and ISAAC (Phillips 
et al. in prep). In the first year of the program, 130 
NIR spectra have been obtained of 23 SNe la. The NIR 
time-series spectroscopy of SN 201 lfe has demonstrated 
the many utilities of NIR spectra, both before and after 
maximum light, for both time-series and single "snap- 
shot" observations. This improved data set will allow us 
to examine the efficacy of these tools and help to further 
decipher the nature of SNe la in the near future. 
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APPENDIX 

PRINCIPAL COMPONENT MODEL LINE PROFILE 

Principal component analysis (PCA) is a statistical technique useful for re ducing the dim ensionality and identifying 
patterns in a set of multi-dimensional data. We use the "bra ket" notation of lSuzukil (|2006T ) here for the formulation of 
PCA. A sample of early, high signal-to-noise ratio NIR spectra from the literature, including those of SN 2011fe, were 
used as PCA data input to build the model of the Mg II A1.0927 /im profile. Each spectrum was shifted in wavelength, 
such that the Mg II line minima coincide at a fiducial wavelength of 1.05 /im, normalized to the same integrated flux 
in the Y band, and interpolated into N fixed wavelength elements. The wavelength interval was chosen to be 10A, 
approximately matching and preserving the medium resolution of the SN 201 lfe spectra at 1.05 fim. The log of the 
flux was used as the PCA input data and is denoted by m). The mean spectrum \fx) was computed and subtracted 
from the data set. A N x N covariance matrix was then computed to determine the strength of the correlation between 
spectral features. 
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The eigenvectors and eigenvalues of the co variance matrix w ere computed using internal IDL routines which utilize 
Householder reductions and the QL method (|Press et al.lll992|) . Each of the N eigenvectors, or principal components, 
denoted here by was ranked by the corresponding eigenvalue, which provides a measure of the amount of data vari- 
ation the principal component describes. The zeroth-ranked principal component points in the direction of maximum 
variance in the AT-dimensiona! data space, and subsequent principal components are orthogonal vectors: 

(Al) 

Each of our input spectra can be expressed as the sum of the eigenvectors: 

N-l 

K) = £>ite>. (A2) 
3=0 

where is the projection of the ith input spectrum on the jth eigenvector. 

Using these eigenvectors, the model spectrum of the Mg II feature was then built as M(p, n), which varies as a 
function of the number of eigenvectors n < N, and projection pj for each eigenvector |£j): 



n-l 

M{n,p) = \ fJ ,) + Y / PMi)- 



(A3) 



